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Abstract. We present results of the performance of the second prototype of the CASTOR quartz-tungsten sampling 
calorimeter, to be installed in the very forward region of the CMS experiment at the LHC. The energy linearity and res- 
olution, as well as the spatial resolution of the prototype to electromagnetic and hadronic showers are studied with E = 
20-200 GeV electrons, E = 20-350 GeV pions, and E = 50, 150 GeV muons from beam tests carried out at CERN/SPS 
in 2004. The responses of the calorimeter using two different types of photodetectors (avalanche photodiodes APDs, 
and photomultiplier tubes PMTs) are compared. 
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1 Introduction 

The CASTOR (Centauro And Strange Object Research) detec- 
tor is a quartz-tungsten sampling calorimeter, which has been 
proposed for the study of the very forward rapidity region in 
heavy ion and proton-proton collisions in the multi-TeV range 
at the LHC Q. Its main physics motivation is to complement 
the nucleus-nucleus physics programme, focused mainly in the 
baryon-free region at midrapidity |2|. CASTOR will be in- 
stalled in the CMS experiment at 14.38 m from the interaction 
point, covering the pseudorapidity range 5.2 < r| < 6.6 and 
will, thus, contribute not only to the heavy ion programme, but 
also to diffractive and low-x physics in pp collisions [3]. The 
results of the beam test and simulation studies with CASTOR 
prototype I |4| prompted us to construct a second prototype 
using quartz plates, avalanche photodiodes (APDs) as well as 
photomultiplier tubes (PMTs), and air-core light-guides with 
inner reflective foil (Dupont polyester film reflector coated with 
AlO and reflection enhancing dielectric layer stack Si02+Ti02). 
In addition, we tested a new semi-octant ((]) = 22.5°) geometry 
of the readout unit in the electromagnetic section. The beam 
tests were carried out in the H2 line at the CERN SPS in 2004 
using beams of electrons, pions and muons. The prototype II 
calorimeter consists of an electromagnetic (EM) and a hadronic 
(HAD) section, built in an octant sector (Fig.Q]). Both calorime- 
ters are constructed with successive layers of tungsten plates 



as absorber and fused silica quartz plates as active medium. 
The EM part (14 cm length) is further divided into two semi- 
octant sectors and is longitudinally segmented into 2 sections, 
so that there are 4 independent readout units in total. The HAD 
part (40 cm length) retains the octant geometry of prototype I 
and is longitudinally segmented into 4 sections. The Cerenkov 
light produced by the passage of relativistic particles through 
the quartz medium is collected in sections along the length of 
the calorimeters and focused by air-core light guides onto the 
photodetector devices, APDs or PMTs. 



2 Technical description 

The CASTOR detector is a Cerenkov-effect based calorimeter 
with tungsten absorber and quartz plates as active material. The 
main advantages of quartz calorimeters are their radiation hard- 
ness, the fast response and the compact detector dimensions [5 ] 
very well adapted for the experimental conditions encountered 
in the very forward region at the LHC. A detailed description of 
the operation principle and, in particular, of the light-guide per- 
formances have been provided in reference (4). In section |2~T1 
we describe the active (quartz) and passive (tungsten) mate- 
rials of the calorimeter considered in this second beam test. 
Section 12.21 discusses the characteristics of the two types of 
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Fig. 1. Picture of the CASTOR prototype II calorimeter before assembling the photodetectors. The semi-octant geometry of the EM section 
(length: 14 cm) and the octant geometry of the HAD section (length: 40 cm) can be seen. 



photodetectors (photomultipliers and avalanche photodiodes) 
tested. 

2.1 Tungsten - Quartz plates 

The calorimeter is constructed from layers of tungsten (W: Xj= 
10.0 cm, Xq= 0.365 cm, density= 18.5 g/cm 3 ) plates as absorber 
and fused silica quartz (Q) plates as active medium (see Fig.[2j. 
For the electromagnetic section, the W-plates have a thickness 
of 3 mm and the Q-plates 1 .5 mm. For the hadronic section, the 
W- and Q-plates have a larger thicknesses of 5 mm and 2 mm, 
respectively. The W/Q-plates are inclined 45° with respect to 
the direction of the impinging particles, in order to maximize 
the Cerenkov light output in the quartz. Each individual combi- 
nation of W/Q-plates is called a sampling unit (SU). The large 
sides of the Q plates were covered with Ty vek paper, to protect 
them from damage by the tungsten plates and also to diffuse 
back the escaping light. The perimeter sides - except the top 
one - were painted with white reflecting paint. The top edge of 
the W plates had just a machined finish. 

In the EM section, each sampling unit (SU) corresponds to 
1.218 X , or4.88xl0~ 2 V Each readout unit (RU) consists of 
11 SUs and is 13.4 Xo, or 0.536 X] deep. The EM section is 
divided in two successive RUs and has a total length of 26.8 
Xq and 1.072A./ lengths. In the hadronic section, a sampling 
unit corresponds to 7.96- 10~ 2 A./. Each readout unit consists of 
10 SUs and is 0.796 A,/ deep. The HAD section has 4 RUs, 
corresponding to 3.186 A/. 

In total, the whole prototype has 4.26 A/. For some runs 
with pions, we inserted an additional inactive absorber of 1 .03 
Xi in front of the calorimeter, in order to make the EM section 
act as a hadronic one, increasing the total depth of the prototype 
to 5.3A/. 

2.2 Photodetectors 

The Cerenkov light emitted by the quartz plates is collected 
and transmitted to photodetector devices through air-core light- 
guides. All light guides of Prototype-II were equipped with 



Dupont [AIO+ Si02+Ti02] reflective foil with the same char- 
acteristics discussed in [4|. As photodetectors we used a ma- 
trix of 4 or 6 Hamamatsu S8148 APDs (developed originally 
for the CMS electromagnetic calorimeter J6|), as well as two 
different types of PMTs. The total area of the APDs was 1 cm 2 
(for 4 APDs) and 1.5 cm 2 (for 6 APDs), see Fig. [3] The pho- 
totubes were positioned only on one side of the EM section of 
the prototype, for comparison with the APDs during the elec- 
tron beam tests. The two types of PMTs used were respectively: 
(i) a Hamamatsu R7899 PMT, and (ii) a radiation-hard multi- 
mesh, small size PMT FEU- 187 from RIE St. Petersburg, with 
cathode area ^2 cm 2 H. 

3 Beam tests 

The beam test of prototype II took place in October 2004 at the 
H2 beam line of the SPS at CERN. Electron, hadron (tit) and 
muon (jU _ ) beams of several energies were used. The energy re- 
sponses (linearity, resolution) of the electromagnetic and hadronic 
calorimeters were obtained with energy scans with: 20-200 GeV 
electrons, 20-350 GeV pions, as well as 50, 150 GeV muons. 
The calorimeter prototype was placed on a platform movable 
with respect to the beam in both horizontal and vertical (x,y) 
directions (see Figure |4j. A telescope of finger scintillator de- 
tectors and wire chambers were installed upstream of the proto- 
type, giving precise information on the position of each particle 
hitting the calorimeter. In this way, we were able to know the 
beam profile and also select particular regions of the beam pro- 
file for the spatial resolution analyses. We note that the typical 
visible transverse sizes of hadronic and electromagnetic show- 
ers in quartz calorimeters are 0(5-10 cm), 0(10 mm) resp. (for 
95% signal containement), i.e. are a factor 3 to 4 times narrower 
than those in "standard" (scintillation) calorimeters 0. 

Figure|5]shows the two semi-octants of the electromagnetic 
(blue) and the octant of the hadronic (red) sections, as seen pro- 
jected onto a plane at 45° with respect to the beam axis. We 
notice that there is no complete overlap of the two sections, 
due to the different sizes of the W/Q-plates available. The hor- 
izontal and vertical numbers correspond to distances along the 
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Fig. 2. Upper photograph of the W/Q-plates of the CASTOR prototype-II showing the EM and HAD sections (lower picture) and the light 
guides (upper picture) in the semi-octant (octant) geometry of the EM (HAD) sections respectively. 



4(6) APD Readout unit* 
with electronics 
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Fig. 3. Assembled APD readout units with 4 and 6 APDs. 
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plate (x — y coordinates) of the points used for the horizontal 
and vertical scans. 



Table Q] lists the (x,y) coordinates of the impact points of 
the horizontal and vertical scans for both electron and hadron 
beams. The location of these points on the 45° projection of the 
semi-octant sectors is shown in Figure [5] The beam profile for 
each point was subdivided into a number of smaller parts, each 
of diameter ~l-2 mm, so that we obtained more impact points 
in total. 



4 Electron beam tests 

Electron beams of energy 20-200 GeV were used to test the 
energy linearity and resolution as well as the position resolution 
of the EM section of the prototype. 

4.1 Energy response 

A typical spectrum measured with 100 GeV electrons incident 
on the EM section of the prototype, equipped with PMTs, is 
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shown in Figure [6] Residual muons in the electron beam are 
also seen as minimum ionizing particle (MIPs) just above the 
pedestal. The energy response of the calorimeter is found to be 
Gaussian for all energies. Figure [7] shows the energy response 
for 20 and 200 GeV electron beams, obtained with 4 and 6 
APDs respectively. 

Energy Linearity: To study the linearity of the EM calorime- 
ter response as a function of electron-beam energy, a central 
point (Fig. [8]) in the two different azimuthal sectors has been 
exposed to beams of various energies. The distributions of sig- 
nal amplitudes, after introducing the cuts on the spatial profile 
of the beam (a circle of radius 2 mm), are in most cases sym- 
metric and well fitted by a Gaussian function. The peak signal 



position, obtained for the three photodetector configurations, is 
plotted as a function of the beam energy in Figure[9] 

For all configurations, the calorimeter response is found to 
be linear in the energy range explored. The average signal am- 
plitude, expressed in units of ADC channels, is satisfactorily 
fitted by the formula: 



ADC = a + bxE (1) 

where the energy E is in GeV. The fitted values of the pa- 
rameters for each configuration are shown in the insets of each 
plot in Fig. [9] 
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Table 1. The (x,y) coordinates (mm) of the impact points of the horizontal and vertical scans for both electron and hadron beams. 
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Energy Resolution: The relative energy resolution of the calorime- Figure[l0]shows the fit to the data with expressions §2^ and 



ter has been studied by plotting the normalized width of the 
Gaussian signal amplitudes, o/E, with respect to the incident 
beam electron energy, £(GeV) and fitting the data points with 
two different functional forms J4]: 



(2) 



o/E = po + pi/VE 



o/E = po® pi/VE © p 2 /E 



(3) 



where the © indicates that the terms are added in quadra- 
ture. In principle, three general terms contribute to the energy 
resolution in calorimeters: 

1 . The constant term, po, related to imperfections of the calorime 
try, signal generation and collection non-uniformity, cali- 
bration errors and fluctuations in the energy leakage, which 
limit the resolution at high energies. 

2. The stochastic or sampling term, p \ , due to intrinsic shower 
photon statistics, characterizes the fluctuations in the signal 
generating process. 

3. The noise term, p2, includes the electronic noise contribu- 
tion from capacitance and dark current which (due to its 
steep l/E dependence) is only important for low energies. 



0. Both parametrizations satisfactorily fit the data. In Table|2] 
we summarize the fit parameters for both parameterizations and 
the three readout configurations. The measured stochastic term 
pi is in the range 36% - 54%. We notice too that the constant 
term po is close to zero for all options. It should be noted that 
though the APDs are very sensitive to both voltage and temper- 
ature changes, there was no stabilization used for this test. 

Spatial Response: The purpose of the area scanning was to 
check the uniformity of the EM calorimeter response to elec- 
trons hitting at different points on the sector area, as well as to 
assess the amount of edge effects and lateral leakage from the 
calorimeter, which could lead to cross-talk between neighbor- 
ing sectors. Figure [8] shows the typical profile of the electron 
beam hitting the left semi-octant of the prototype. The width 
of the EM shower and the percentage of the containment close 
to the edge were estimated by varying the horizontal and ver- 
tical hit positions of the incident beam according to the (x,y) 
coordinates shown in Fig. [8] and listed in Table Q] 

The results of the horizontal-scan analysis are shown in 
Figure [TT] for the 4 APDs readout configuration. Figure \TT\ a) 
shows the response of the two adjacent (left-right) EM semi- 
octants as the beam impact point moves across the front face 
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Electrons @ 100 GeV with PMT's 
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Fig. 6. Energy response of the EM calorimeter equipped with PMTs to 100 GeV electrons (and residual beam muons). 



Table 2. Energy resolution parameters of the EM calorimeter prototype as obtained from the measured electron beam energy resolution and 
Eqs. (H, (H. 
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of the calorimeter. The sigmoid nature of each response curve 
is evident. In Figure fTIT b). the x— derivative of the response is 
calculated, giving the width of the electromagnetic shower. We 
observe that one standard deviation amounts to 1.7 mm. 



The vertical-scan covered the entire height of the semi- 
octant EM sector, with impact points shown in Figure and 
listed in Table Q] The results of this scan are shown in Fig- 
ure Q~2] We notice the abrupt fall at the lower end of the sector 
past the point "A" and the more gradual fall at the upper end, 
the later due to the shower particles directly hitting the light 
guide. 



4.2 Pion beam tests 

Pions of energy 20-350 GeV were used for the study of the 
hadronic energy and position responses of the CASTOR proto- 
type II. In order to increase the interaction depth of the calorime- 
ter, an inactive absorber of 1.03A./ was inserted in front of the 
EM calorimeter, increasing the total depth to 5.3A,/. This had 
also as a result to make the two first (EM) RUs effectively act, 
in depth, as part of the hadronic section. 

Energy Response: Typical spectra, obtained with 200 GeV 
pions incident on the prototype, are shown in Figure [T3l where 
the distribution of the total energy measured in both (EM and 
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4 APD's, electrons @ 20 GeV 4 APD's, electrons @ 200 GeV 
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Fig. 7. Energy response of the EM calorimeter to electron beams of 20 and 200 GeV obtained with 4 APDs (upper plots) and 6 APDs (bottom 
plots). 



HAD) parts of the calorimeter is plotted. During the different 
tests, the electromagnetic sections were equipped with 4 or 6 
APDs and the hadronic ones had 4 APDs in its readout units 
for all runs. The total depth of the prototype (5.3A./) was not 
enough to contain the showers produced by the pion beams. 
We see that there is a long tail at high energies indicating the 
leakage of energy from the back of the calorimeter. However, 
no quantitative measurements of the leakage fraction were done 
at this stage. The peak of the total pion energy measured by 
the prototype was fitted with a Gaussian and a Landau curve. 
The fitting ranges correspond roughly to 1- (2-)o around the 
peak for the Gaussian (Landau) distributions. We observe that 
the Landau parametrization fits the distribution better than the 
Gaussian one. 

The energy response (position and width of the pion peak) 
was obtained by fitting both Gaussian and Landau curves to the 
spectrum measured for all beam energies. The corresponding 
hadronic energy linearity and resolution were thus obtained. 

Energy Linearity: Figure [14] shows the linearity of the CAS- 
TOR prototype to incident pions as obtained by measuring the 
total energy deposited in the calorimeter sections and corre- 
lating the position of the pion peak with each corresponding 
beam energy. At higher energies, the Landau fit gives higher 
response, as expected, and an overall smaller statistical error. 



Energy Resolution: The relative energy resolution of the calo- 
rimeter has been studied by fitting the normalized width of the 
fitted signal amplitudes (peaks in Fig.QjJl, o/E, with respect to 
the incident pion beam energy, £XGeV), with the two functional 
forms d2J an d ®- Figure[T5]shows the obtained energy resolu- 
tion of the prototype for pions of energy up to 350 GeV with 
6 (left) and 4 (right) APDs in the EM part of the calorimeter. 
The blue points and line in Figure[l5]show the resolution when 
the pion energy distribution is fitted by a Gaussian curve. The 
red ones, when the distribution is fitted by the Landau expres- 
sion. We observe that the resolution is much better when the 
Landau fit is employed. It should be noted that the length of 
the tested calorimeter is only 4.26 interaction lengths (almost a 
factor 3 smaller than the planned length of the final CASTOR 
calorimeter) and thus there is considerable energy leakage at 
the end even at low pion energies. This does not permit an ac- 
curate estimation of the hadronic resolution. 

Spatial Response: Figure[l6]shows the pion beam profile hit- 
ting the left semi-octant region of the prototype. We observe 
that the hadron beam is much more focused than the electron 
beam (see profile in Fig. [8]). The spatial response of the proto- 
type calorimeter to pions was obtained from the two EM semi- 
octant sectors, by moving the beam along the x— direction. The 
1 .03 Xj inactive absorber was positioned in front of the calorime- 
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Beam Profile, electrons @ 200 GeV, 4 APD's 
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Fig. 8. Profile of 200 GeV electron impinging on the left semi-octant of the calorimeter, as measured by the scintillator-wire-chamber telescope 
upstream of the prototype. 




Fig. 9. Energy response linearity (signal peak-position versus beam energy) of the EM section, obtained with different photodetectors: 4 APDs 
(left), 6 APDs (center), and PMTs (right). 
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Resolution of CASTOR Proto II: c" , 4 APD's (3 parameter fit) 
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Fig. 10. Energy resolution (signal peak width versus beam energy) of the prototype EM section, obtained with the three readout configurations 
considered: 4 APDs (left), 6 APDs (center), and PMTs (right): 3-parameters fit Eq. l|3} (top); 2-parameters fit Eq. {2]l (bottom). 
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Fig. 11. (a) Response of the left and right semi-octant sectors of the EM section as the beam scans the front face of the calorimeter, (b) The 
derivative of the response with respect to x, indicating the width of the EM shower. 
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Y Area Scan, e @100 GeV 
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Fig. 12. Spatial scan along the y-direction of one EM sector. The impact points are those listed in Table Q] and shown in Fig. [5] [Note that, at 
variance with the beam-coordinate-system used in Fig. [5] we use here a calorimeter-coordinate-system, based on the position of the calorimeter 
on the moving platform]. 
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Fig. 13. Total energy spectra (ADC channel counts) measured in the prototype-II for the pion beam of 200 GeV and 6 APDs in the EM section. 
The pion peak is fitted to a Landau (top plot) and Gaussian (bottom plot) curve with fit parameters reported in the inset. The peak to the left is 
the pedestal. 
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Linearity of Castor Proto II: Pions with 4 APD's (EM + HA) 
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Linearity of Castor Proto II: Pions with 6 APD's (EM + HA) 
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Fig. 14. Energy response linearity of the prototype calorimeter to pions of several energies, fitted to Gaussian (blue) and Landau (red) 
parametrizations. The top (bottom) plot is obtained with 4 (6) APDs readout in the EM section. 



ter. The beam profile for each point was subdivided into a num- 
ber of parts, each of diameter ~5 mm, so that we obtained more 
impact points in total. FigurefTTlshows the a— scan for pions of 
300 GeV energy on the left and the derivative of this response 
with respect to x on the right. The pion beam width has Chad 
= 6.4 mm, considerably larger than the corresponding electro- 
magnetic one (Gem = 1-7 mm, see Fig. fTTb. as expected. 



5 Muon beam tests 

Muon energy spectra at 50 and 100 GeV were measured with 
the electromagnetic sector, using the PMT readout configura- 
tion. Figure [l8]shows the muon peak measured for the 50 GeV 
beam well separated from the pedestal at zero counts. The line- 
shape has been obtained with two different PMTs: Hamamatsu 
R7899 (Fig. [HJi), and RIE FEU187 <H3>). In Figure [18};, the 
sum of both EM readout units is shown. 

From Figure [18] we find that the Hamamatsu R7899 PMT 
performs much better than the RIE FEU 187 one, in identify- 
ing the muon signal above the pedestal. A disadvantage of the 
R7899 PMT for this application is its large length, which pro- 
hibits its use, even in the semi-octant geometry. 



6 Monte Carlo simulation of prototype II 



Fig. [19] shows the GEANT4 [7 1 geometry of prototype II as 
implemented in the CMS software (OSCAR 6.3.5). The geom- 
etry of the electromagnetic section described in the simulations 
(XML-format) matches exactly that of the tested calorimeter. 

We run simulations for 500 electron events with 7 differ- 
ent energies in the range E = 20 -250 GeV and studied the 
corresponding response in terms of the number of photoelec- 
trons produced. Figure [20] shows the simulated energy (a) lin- 
earity and (b) resolution of the prototype obtained assuming an 
overall efficiency (light transmission x quantum efficiency) of 
about 65% for the APDs |4|. The linearity of the energy re- 
sponse is consistent with the experimental data (Fig.[9jl, but the 
energy resolution is 2-3 times better than the beam test results 
(Fig.EB. 

Figure [22] shows the jt-spatial response of the electromag- 
netic shower simulated in GEANT4. In the MC simulation, the 
electron beam has a radius of 1.5 mm, similar to the cut im- 
posed in the analysis of the experimental data. The sigmoid 
curve is seen in Fig. l22k and its ^-derivative in Fig. [22b, from 
which we obtain the width of 1 .56 mm. which is close to what 
one observes in the real data (Fig.fTTb) 
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Resolution of Castor Proto II: Pions with 4 APD's 
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Resolution of Castor Proto II: Pions with 6 APD's 
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Fig. 15. Energy resolution, o/E, of CASTOR prototype II to pion beams of several energies obtained with 4 (top) and 6 (bottom) APDs readout 
in the EM section. The different fit parameters shown in the inset are obtained with Eq. (f3j and Eq. (O when the widths o of the pion peaks are 
fitted to a Gaussian or Landau distribution. 



7 Summary 

We have presented a detailed performance study of the ener- 
getic and spatial responses of a second prototype of the CAS- 
TOR quartz-tungsten calorimeter of the CMS experiment. The 
results have been obtained from beam tests at CERN-SPS with 
high-energy electrons (20-200 GeV), pions (20-350 GeV) and 
muons (50, 150 GeV) and two different types of photodetectors 
(APDs and PMTs) for the EM section of the calorimeter. The 
main conclusions of this study can be summarized as follows: 

1 . EM Section: The semi-octant geometry has an efficient light- 
collection with 4 or 6 APDs. Due to the small height of the 
light-guide, a PMT readout can also be used, provided it is 
of small size. This has the advantage of higher gain (over 
the APD configurations), enabling the clear identification 
of the muon peak above the pedestal. 

2. HAD Section: The octant geometry has an efficient light- 
collection for the hadronic section. However, the large height 
of the associated light-guides precludes this configuration 
in the limited space available for the CASTOR calorimeter 
in the very forward region of the CMS experiment. 

On the basis of physics concerns for both pp and heavy- 
ion interactions, the semi-octant geometry (which would cor- 



respond to 16 sectors covering full (])) is, therefore, preferred. 
For this geometry, two reading-device options provide the de- 
sired performances: (i) 6 Hamamatsu-S8148 APDs per readout 
unit, and (ii) a small-size PMT, such as the RIE FEU- 187. Both 
photodetectors should be tested/adapted for the radiation-harsh 
conditions of the CASTOR calorimeter (~ 10-100MGy accu- 
mulated, to be compared e.g. to the ^3kGy expected for the 
CMS ECAL APDs JU). The relative merits and difficulties of 
each option will be further studied in detail before a final deci- 
sion is reached. 
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Fig. 16. Profile of the 300 GeV pion beam impinging on the left semi-octant region of the calorimeter. 



X Area Scan, pions @ 300 GeV 



Derivative Curve X-Scan: pions @ 300 GeV 
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Fig. 17. x— scan along the face of the prototype for 300 GeV pions (left plot). The derivative of the sigmoid curve, giving the width of the 
hadronic shower distribution (right plot). 
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Fig. 18. Energy spectra measured in the EM section of prototype II with a muon beam of 50 GeV energy and using two different PMTs: (a) 
Hamamatsu R7899, (b) RIE FEU187, and (c) the sum of both. 



MC Simulated Prototype Geometry 




Fig. 19. Layout of the simulated geometry of the CASTOR prototype as implemented in GEANT4 (OSCAR 6.3.5). 



Linearity of Castor Prototype II 



Resolution of CASTOR Prototype II 
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Fig. 20. Simulated energy response in terms of photo-electrons generated in the EM sections of the CASTOR prototype: (a) linearity, (b) 
resolution. 
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Resolution of CASTOR Proto II, electrons 
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Fig. 21. Comparison of the experimental resolution for the three light readout configurations considered (4,6 APDs and PMTs) and the MC 
simulated one. 



X Area Scan, electrons @ 100 GeV 



Derivative Curve X-Scan: electrons @ 100 GeV 
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Fig. 22. (a) Simulated x— profile of the electromagnetic shower, (b) Derivative of the simulated response with respect to x, indicating the width 
of the EM shower. 
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